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Gamma-ray polarimetry is a new and prospective tool for studying various extreme high-energy celestial 
objects and is of great significance for the development of astrophysics. With the rapid development of mi- 
crosatellite technology, the advantages in space exploration are becoming increasingly apparent. Therefore, in 
this paper, we conducted a simulation study on a soft gamma-ray polarimeter for a microsatellite in space. Here, 
we proposed a unique design structure for the polarimeter based on the microsatellite design concept and the 
Compton scattering principle. And then, the detailed Monte Carlo simulations using mono-energetic gamma-ray 
linear polarization sources and the Crab-like sources in the energy range of 0.1-10 MeV with full consideration 
of the orbital background were performed. The simulation results show that the polarimeter can exhibit ex- 
cellent polarization detection performance. The modulation factor is 0.80-+0.01, and the polarization angles 
are accurate within an error of 0.2° at 200 keV for on-axis incidence. For the Crab-like sources for on-axis 
incidence, the polarization degrees are consistent with the set values within the error tolerance, the modulation 
factor is 0.76+0.01, and the minimum detectable polarization reaches 2.4% at 30 for an observation time of 
10° seconds. In addition, the polarimeter has recoil electron tracking, imaging, and powerful background sup- 
pression at a large field of view (~27 sr). The polarimeter designed can meet the requirements of a space-soft 


gamma-ray polarization detector very well and has a bright research prospect. 


Keywords: Soft Gamma-ray, Polarization, Compton scattering, Microsatellite, Monte Carlo simulation 


I. INTRODUCTION 


Gamma-ray polarization measurement in gamma-ray as- 
tronomy is broadly considered to be a new and powerful 
diagnostic tool for some critical open questions or doubts 
about the most extreme high-energy sources: gamma-ray 
bursts (GRBs), pulsars, active galactic nuclei (AGNs), bi- 
nary black holes (BBHs), etc. [1, 2], and these questions or 
doubts might not be resolved or explained through the tim- 
ing, energy, and direction of gamma rays. For GRBs studies, 
polarization measurements can illuminate the nature of the 
central engines that produce ultra-relativistic jets in GRBs, 
as well as the physical properties, the radiation mechanisms, 
and the energy dissipation points of these jets, and can also 
contribute to the constraints of the theoretical models of the 
origin of GRBs [1, 3-6]. Detecting gamma-ray polarization 
emitted by pulsars can help estimate the magnetic field struc- 
ture around compact objects, to understand the mechanism 
of gamma-ray emission (curvature radiation or synchrotron 
radiation), and to speculate on particle acceleration and pair- 
cascading processes in the magnetosphere of pulsars [1, 7]. 
In particular, accreting black hole (BH) systems, including 
BBHs and AGNs, are considered to radiate linearly polar- 
ized X-rays and gamma rays owing to scattering processes 
in their accretion disks, and therefore, the measurement of 
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these polarization features will allow us to identify the ge- 
ometry of the corona [2, 8]. Many more examples of space 
gamma-ray polarization as a unique tool to study high-energy 
astrophysics can be found in Ref. [9-13]. In conclusion, 
gamma-ray polarization is uniquely helpful and valuable for 
astrophysical studies as an up-and-coming tool, which makes 
space gamma-ray polarization measurements extremely at- 
tractive. 


Under the tremendous scientific lure of space gamma- 
ray polarization, many astrophysicists from various coun- 
tries have joined the team to detect space gamma-ray po- 
larization and have made efforts in polarization detection 
by launching satellites and flying high-altitude balloons. A 
few typical polarization-related satellite experiments are an 
X-ray polarimeter onboard the eighth Orbiting Solar Obser- 
vatory (OSO-8) mission [14], the Ramaty high energy so- 
lar spectroscopic imager (RHESSI) mission [15], IBIS and 
SPI onboard International Gamma-ray astrophysics labora- 
tory (INTEGRAL) [12, 16], GAP onboard IKAROS space- 
craft [17], POLAR placed on the Chinese space station TG- 
2 [18], and cadmium zinc telluride imager (CZTI) onboard 
ASTROSAT satellite [19]. In addition to the satellite ex- 
periments mentioned above, there are several representative 
balloon experiments: the Polarimeter for High ENErgy X- 
rays (PHENEX) [20], the Polarised Gamma Ray Observer 
(PoGO) [21], PoGO+ [22], the Gamma RAy Polarimeter Ex- 
periment (GRAPE) [23], the Gamma-Ray Astro-Imager with 
Nuclear Emulsion (GRAINE) experiment [24]. Although the 
number of polarization detection experiments available is rel- 
atively large (as listed above), most are focused on the X- 
ray band, especially the hard X-ray band (tens to hundreds 
of keV), with a few detecting high-energy (>~10 MeV) 
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gamma-ray polarization, and an extreme lack for the soft 
gamma-ray band (~0.1-10 MeV), specifically above 1 MeV. 


Generally, the most efficient way to detect the polariza- 
tion of soft gamma rays is based on the Compton scatter- 
ing principle. This is because the radiation physics pro- 
cesses occurring in soft gamma rays with matter are domi- 
nated by Compton scattering, similar to detecting low- and 
high-energy gamma-ray polarization dominated by the pho- 
toelectric effect and the electron-positron pair production, re- 
spectively. Here, it is worth noting that the Medium Energy 
Gamma-Ray Astronomy (MEGA) telescope [25], based on 
the Compton scattering principle, has been shown in ground 
tests to be capable of detecting soft gamma-ray polarization. 
Unfortunately, the MEGA project was not completed for var- 
ious reasons. In addition, the Compton Spectrometer and Im- 
ager (COSI) [26], also based on the Compton scattering prin- 
ciple, is equally sensitive to gamma-ray polarization in this 
energy band. However, COSI is only an end-of-flight balloon 
experiment and faces many challenges in gamma-ray polar- 
ization detection as a balloon payload. The most important of 
these challenges are atmospheric absorption and scattering, 
and limited exposure. This significantly reduces the ability of 
the detector to detect soft gamma-ray polarization in space. 
The best solution to the above problem is to send the detector 
into space, namely by launching a detection satellite. Given 
the above analysis, there is a great need to propose and im- 
plement a satellite project on soft gamma-ray polarimetry in 
space to occupy this almost empty energy region, which is of 
great scientific importance for astrophysical research. 


There is no doubt that the best way to measure gamma-ray 
polarization in the universe is to launch detection satellites. 
What needs to be highly concerned in recent years is that 
with the development of science and technology, the develop- 
ment of modern small satellites has become more and more 
rapid and has significant advantages: small size, lightweight, 
high technology, short development cycle, low cost, can be 
standardized stars and modular design technology, can be 
mass production and storage in the flow line and easy to 
launch [27, 28]. However, compared with modern small satel- 
lites, traditional exploration satellites have apparent disadvan- 
tages such as large size and weight, complex technology, long 
development period, high cost, high risk, and difficulty in 
achieving [27, 28]. Given this, modern small satellites can 
give a new perspective to cosmic gamma-ray polarimetry. 


To achieve a high-precision and all-sky survey of the po- 
larization of soft gamma rays in space, the idea of using a 
constellation composed of multiple microsatellites is promis- 
ing, like the ”GRID mission” [29]. Apparently, excellent po- 
larization detection performance is essential for each of the 
microsatellites in the constellation. Therefore, in this paper, a 
detailed study was carried out for a microsatellite (polarime- 
ter): a novel detector design based on the principle of Comp- 
ton scattering was proposed to detect the linear polarization of 
soft gamma rays, a detector mass model used for simulation 
was constructed, and a detailed simulation experiment using 
Monte Carlo methods was implemented to verify its perfor- 
mance. The simulation results showed that the polarimeter 
designed has excellent performance. Overall, this work vali- 
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dates the rationality of our proposed detector design scheme, 
lays a solid foundation for future polarimeter development, 
and offers a meaningful reference for other researchers to de- 
sign and develop detection satellites. 


II. PRINCIPLES AND METHODS 


A. Principle of polarization measurement 


In the soft gamma-ray domain, the Compton scattering ef- 
fect becomes a dominating process in interactions of photons 
with matter. Compton scattering preserves the polarization 
information of linearly polarized photons to a certain degree. 
Expressly, when linearly polarized photons and matter un- 
dergo Compton scattering, the azimuthal angle distribution of 
the scattered photon is related to the polarization degree and 
direction of the incident photon. Therefore, the soft gamma- 
ray polarimeter in this paper is based on this principle to de- 
tect the gamma-ray polarization signals in the universe. 

If the initial gamma ray is linearly polarized, then its 
Compton scattering process can be graphically demonstrated 
in Fig. 1. Additionally, the Compton scattering differen- 
tial cross-section of the linearly polarized photon can be ex- 
pressed by the Klein-Nishina formula as [18] 


do ri (E ?2/E E 3 ‘ 
= 2 s 8 
dQ 5) (=) E T E sin^ @ cos* 7) 
2 IN 2 
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where ro is the classical electron radius, Æ is the incident or 
initial photon energy, F’ is the scattered photon energy, 0 is 
the Compton scattering angle, and 77 is the angle between the 
scattering direction of the scattered photon and the polariza- 
tion direction of the incident photon (i.e., the azimuthal angle) 
shown in Fig. 1. E’/EF in Eq. (1) can be represented by the 
Compton equation [1] 


(2) 


where mec? is the electron rest mass energy. It can be seen 
from Eq. (1) that linearly polarized photons tend to scatter 
perpendicularly to the incident polarization vector (minimiz- 
ing the term 2sin? 4 cos? n). In addition, when both E’ /E 
and 0 are constants, the differential scattering cross-section 
of polarized photons will obey the cos[2(7 + 2/2)] distribu- 
tion, namely, the variation of the number of scattered photons 
with the azimuthal angle ņ obeys the cosine distribution. In 
practice, the incident photon energy E and the scattering an- 
gle 0 are often non-constant but take on a range of values, 
and what we actually measure is the averaging effect, but this 
does not affect the fact that the averaged azimuth follows a 
cosine distribution. 
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Fig. 1. (Color online) The schematic of Compton scattering of a 
linearly polarized photon. 


In most cases, since 7 cannot be measured directly, the po- 
189 larization information of incident linearly polarized photons 
) 160 cannot be obtained using the method of measuring 7. Given 
) 161 the above practical situation, the angle ọ is introduced as the 
) 162 Compton scattering azimuth, which is the angle between the 
| 13 polarized scattered photon plane and the X-axis, as shown in 

64 Fig. 1. Polarization signatures of a source of incident linearly 
| 165 polarized photons are reflected in the distribution of the az- 
imuthal angle ọ that can be described by the following func- 
tion 


f(y) =A{1 + woos [2 (p — po) + 7]} , (3) 


which can be easily deduced from Eq. (1). Here, yo rep- 
resents the polarization angle of the incident photon or the 
direction of the original polarization vector (see Fig. 1), A 
is the offset of the distribution of the azimuthal scatter angle 
173 (see Fig. 2), and ys is a significant parameter called modula- 
174 tion factor which can describe the polarization response of a 
polarimeter. Equation (3) can be shown graphically in Fig. 2, 
176 and the curve in Fig. 2 is called the modulation curve. The 
177 modulation factor u can be expressed by the equation 
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178 H= 
177 Which can be derived from Eqs. (1), (2), and (3), where Fmax, 
130 Fmin, A, and B are represented in Fig. 2. For a fully (100%) 
181 linearly polarized photon beam, the modulation factor is ex- 
182 pressed as 4100, While for a photon beam with an unknown 
183 polarization degree, the modulation factor is u. The polar- 
184 ization degree of the incident photons can be obtained by the 
185 equation 


(5) 
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187 Where P is a positive value between 0 and 1. 
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Fig. 2. (Color online) The distribution of the azimuthal scatter angle 
y. 


188 Figure 3 illustrates the relationship (the functional relation- 
189 Ship expressed in Eq. (4)) among the modulation factor 4100, 
19 a key performance parameter of a polarimeter, the Compton 
191 Scatter angle 0 and the incident photon energy E. It can be 
192 intuitively seen from Fig. 3 that the modulation of the az- 
193 imuthal distribution is most significant at lower energies and 
194 medium scattering angles. For photons with extremely large 
195 (0 180°) or small (0 0°) scattering angles, they carry little 
196 polarization information. Additionally, we can also see that 
197 the modulation factor 410o decreases significantly with the 
198 increase of photons energy. 


Modulation factor Hoo 


Fig. 3. (Color online) The modulation factor 4100 as a function of 
the Compton scatter angle 0 and the incident photon energy FE. 


199 Besides the modulation factor, the minimum detectable po- 
200 larization (MDP) (also known as polarization sensitivity, i.e., 
201 the detection limit of the degree of polarization [30]) is also an 
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essential performance parameter to describe the performance 
of a polarimeter. The MDP is used to judge the polarization 
detection capability of a polarimeter, which can be calculated 


as follows [1] 
IS+B 
T ? 


where no is related to the expected confidence level of the 
detection (e.g., 2, = 3), and S and B are the count rates of 
the source and background (after all event selection cuts are 
applied) in the observation time T, respectively. As can be 
seen from Eq. (6), the MDP is related to five parameters. In 
general, no, H100, and T are all constants once the confidence 
level, the polarimeter, and the observation time have been de- 
termined. Then the MDP is only affected by S and B, and the 
MDP gets better with increasing © and worse with increasing 
B. Considering that the angular resolution will be used as the 
key event selection condition when calculating S and B, S 
and B are then related to the angular resolution of the detec- 
tor. In fact, the excellence of angular resolution has almost no 
effect on S, while it has a significant effect on B. When the 
angular resolution is excellent, B becomes smaller, resulting 
in a better MDP for the polarimeter. In contrast, when the an- 
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223 gular resolution is poor, B becomes larger, leading to a worse 
224 MDP for the polarimeter. From the above analysis, it can be 
seen that the angular resolution significantly affects the polar- 
ization sensitivity of a polarimeter, which will be a valuable 
guide to the design of a polarimeter. 
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228 B. Polarimeter design 


29 In this paper, the polarimeter was designed to detect the 
polarization information of linearly polarized gamma rays in 
the energy range of 0.1-10 MeV based on the Compton scat- 
tering principle. The structure design of our proposed po- 
larimeter is shown in Fig. 4, where the whole model contains 
only active materials. As a whole, the polarimeter is mainly 
made up of three detection sub-systems: a silicon converter 
(blue) located in the upper center of the polarimeter, a CsI 
absorber (red) surrounding the converter on five sides (ex- 
cept the top side), and an organic plastic scintillator antico- 
incidence shield (ACS) (green) which envelops the two sub- 
detectors mentioned above. The whole detector containing 
only sensitive material has a size of 26x26x20 cm? and a 
mass of ~15 kg. The design details and description of each 
sub-detector of the polarimeter are given in the following. 
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TABLE 1: Design parameters of each sub-detector in the polarimeter. 


Photon Number of 
Design Sensitive Total mass Cell size © Numberof absorption . strips on 
Sub-detectors ; 3 a spacing 
shape materials (kg) (cm”) cells probability (cm) a cell 
@ 1MeV (%) and strip pitch 
Converter a Silicon ~0.2 10x10x0.1 10 ~15 0.5 50,2mm 
Pixel-type Bar:1x1x6 Bar:256 Bar:~80 
Absorber array cal LLS  Cube:2x2x2 Cube:128  Cube:~42 j 
ACS Hollow shell Organic plastic 33 26x 26x20, 1 - i i 


scintillator 


Converter The converter plays a vital role in the po- 
larimeter. It performs two main tasks: first, the first Compton 
interaction takes place in the converter, and second, the con- 
verter records the deposited energy and the interaction po- 
sition of all generated particles. To accomplish the above 
tasks, the converter is required to enable photons to have 
a high Compton scattering probability in it, to be able to 
stop recoil electrons, and to have better energy and posi- 
tion resolution. After comprehensive consideration, it was 
agreed that a multi-layer double-sided silicon strip detector 
was the best converter design choice. Silicon, a low-Z ma- 
terial, has a higher Compton cross-section than medium-Z 
and high-Z materials such as Ge and CdZnTe in the 0.1 to 
10 MeV energy range. The double-sided silicon strip detec- 
tor has a low threshold and guarantees excellent energy and 
two-dimensional (i.e., X and Y direction) position resolution. 
Furthermore, multiple thin-layer configurations can enlarge 
Compton scattering cross section, provide longitudinal posi- 
tion information, and track and absorb recoil electron energy 


thickness: 1cm 


Anticoincidence detector 


Absorber 


Fig. 4. (Color online) The wireframe model (Left) and the mass 
model (Right) of the polarimeter, which consists of the outer organic 
plastic scintillator anticoincidence shield (green), a silicon converter 
(blue) located in the upper middle, and a CsI absorber (red) sur- 
rounding the five sides of the converter. 


24 as much as possible. The overall geometric model of the con- 
25 verter and its physical location in the overall system can be 


seen in blue in Fig. 4, and the detailed design parameters are 
given in Table 1. 


Absorber The absorber (also known as a calorimeter) 
is required to stop scattered Compton photons and measure 
their energy and position information. Additionally, since 
the Compton scattering polarization signature is particularly 
significant at larger scattering angles and lower energies, the 
absorber must be able to measure large-angle scattered pho- 
tons. Finally, the absorber is expected to act as a barrier to 
reduce the radiation from the space orbit environment to the 
converter and improve the background rejection of the entire 
detector. Therefore, an absorber that meets the requirements 
needs to use high-Z materials as the detection medium, has 
good energy and position resolution, has a large acceptance 
of scattered photons, and can surround the converter as much 
as possible. The CsI scintillator is the better choice as the 
sensitive material of the absorber because of its high density, 
large atomic number, high light yield, good mechanical prop- 
erties, not easy to deliquescence, high detection efficiency, 
ease of processing into small pixels of various shapes, avail- 
ability in large quantities, and reasonable price [31]. A pixel- 
type absorber composed of CsI crystals was the solution we 
finally chose. Two types of crystal cells were used: cubic 
crystals and bar crystals. The overall geometric configuration 
of the absorber and its physical position in the overall system 
is given by the red part in Fig. 4, and some design param- 
eters can be found in Table 1. Considering the power con- 
sumption and space constraints of a compact microsatellite, 
silicon photomultipliers (SiPMs) were planned to be used as 
photoelectric converters for CsI crystals instead of traditional 
photomultiplier tubes (PMTs) owing to their attractive abili- 
ties, such as their super miniature size, low power consump- 
tion, low weight, small size, fast time response, large self-gain 
(10°-10°), high signal-to-noise ratio and insensitivity to mag- 
netic fields [29, 32]. Each small cubic crystal was expected 
to couple the photoelectric converter only on the side facing 
away from the silicon converter to minimize the passive ma- 
terial between the converter and absorber. For bar crystals, a 
dual-ended readout scheme will be used. To verify the feasi- 
305 bility of the scheme, we have done a detailed test study for the 
CsI detection cell in the laboratory, the details of which can be 
37 read in Refs. [31, 33]. The experimental results have shown 
that the CsI detection cell exhibits good performance, the en- 
ergy resolution is close to 5% (full width at half maximum, 
FWHM), and the longitudinal position resolution is approx- 
imately 5 mm (FWHM) for 662 keV gamma-ray emitted by 
the 137Cs source. The method of the dual-ended readout can 
not only ensure better energy resolution but also give the posi- 
tion information along the crystal bar direction (Z-direction). 


ACS The ACS of the polarimeter is mainly to veto back- 
ground events induced primarily by charged particles (e.g., 
protons, alphas, electrons, and positrons) from the orbital en- 
vironment in space (described in Sec. IID). A shell-shaped 
hollow plastic scintillator with a thickness of 1 cm was used 
as an anticoincidence shield to completely cover the converter 
and the absorber. The green part in Fig. 4 shows the geometric 
configuration of the ACS and its physical location in the over- 
all system. Meanwhile, some design parameters are given in 
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324 Table 1. 


C. Simulation and analysis tools 


3 The Medium Energy Gamma-ray Astronomy library (ME- 
327 GAlib) [34] is an open-source Monte Carlo simulation and 
ses data analysis package wholly written in C++ and based on 
32 ROOT [35] and Geant4 [36]. It was explicitly designed for 
330 gamma-ray detectors in the low-to-medium energy region. It 
331 could be used for the design of the detector geometry and 
the simulation of the interaction process of gamma rays and 
333 Other particles with matter, as well as data analysis. The 
334 reliability of the MEGAIib package for the simulation and 
data analysis of low- and medium-energy gamma-ray detec- 
tors has been recognized by researchers in this specialized 
field. To date, MEGAIib has been successfully applied to var- 
ious hard X-ray/gamma-ray telescope projects and studies in 
space and on the ground, such as MEGA [37], COSI [26], 
AMEGO [38], COMPTEL [39], ACT [40], TIGRE [41], e- 
ASTROGAM [7], a combined Compton and coded-aperture 
telescope for medium-energy gamma-ray astrophysics [42], 
and many more. 

The simulation of a detector based on MEGAIib first re- 
quires the creation of a realistic geometry by using the Ge- 
omega package contained in MEGAIib [34, 37, 41], as shown 
in Fig. 4. The geometry includes the shape, size, location, 
material, and properties of the surrounding environment for 
each volume that makes up the detector. The cosmic simu- 
lator, Cosima, integrated with the MEGAIib package, is used 
to perform Monte Carlo simulations [34, 37, 41]. Cosima can 
combine Geant4 and the specified source to simulate particle 
transport and interaction with geometric materials and then 
generate an output file that stores the simulated interaction 
information. The simulation data analysis is performed using 
Revan and Mimrec contained in MEGAIib [34, 37, 41, 43]. 
Then, we can get the simulation results. In addition, ME- 
GAlib needs to use the background data computed and gener- 
ated by the LEOBackground software package [44], written 
entirely in Python, to simulate the low-earth orbit (LEO) envi- 
ronment. The above is only a brief overview of the functions 
of the main sub-packages of MEGAIib; a detailed description 
and usage can be found in Ref. [34, 37, 41, 43, 44]. 


332 


D. Simulation method 


Currently, simulation experiments are the primary method 
we use to study the performance of our designed polarimeter. 
In this section, the simulation method of the polarimeter is il- 
38 lustrated, including an overview of the simulation experiment 
39 flow, the detailed configuration of the performance param- 
370 eters of each sub-detector, the configuration of the particle 
371 sources (gamma-ray sources and background sources), and a 
372 brief data analysis process. 

33 Figure 5 shows a brief flow of our simulation experiments 
374 using the MEGAIlib package. First, a realistic geometric 
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model of the detector was constructed, while reasonable per- 
formance parameters were set for the detector in order to 
adapt the simulation to reality. In addition, gamma rays and 
background sources were similarly configured. It should be 
mentioned that the configuration of particle sources also in- 
volved setting information such as the physics list, data out- 
put formats, data storage files, simulation stop conditions, etc. 
Next, Geant4 would be called by MEGAIib to complete the 
Monte Carlo simulation. The simulated data were then recon- 
structed (i.e., Compton event reconstruction) using the proven 
event reconstruction algorithm in MEGAIib. Finally, the sim- 
ulation results were obtained after a high-level analysis of the 
reconstructed events using MEGAIib. A brief description of 
the simulation process has been given above, and a more de- 
tailed elaboration will be shown in the following. 
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Fig. 5. (Color online) A brief flow of simulation experiments using 
MEGAIib. 


As a realistic geometric model of the polarimeter con- 
structed using the MEGAIib package has been shown in Fig. 4 
and described in detail in Sec. IIB, it will not be repeated 
here. However, it is essential to elaborate on the performance 
parameter settings of the detector. To make the simulation 
results reflect the actual situation of the detector as much 
as possible, the setting of simulation parameters is crucial, 
which means that reasonable parameters can make the per- 
formance of the detector obtained by the simulation exper- 
iment closer to that of the real detector and more reliable. 
The double-sided silicon strip detectors that made up the con- 
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verter were assumed to have a uniform energy resolution of 
10 keV FWHM as well as a noise threshold of 15 keV and 
a trigger threshold of 30 keV [37]. The position resolution 
of the converter was determined by the number of strips and 
the thickness of each silicon wafer (described in Sec. I B). 
For the CsI absorber, energy resolutions of 15% FWHM at 
100 keV, 9% FWHM at 350 keV, 6.5% FWHM at 511 keV, 
5% FWHM at 662 keV, 3.5% FWHM at 1000 keV and 2.7% 
FWHM at 5000 keV, a noise threshold of 30 keV, and a trigger 
threshold of 50 keV were set [31, 33, 37]. The depth resolu- 
tion (Z-direction) of the 6 cm CsI crystal was assumed to be 
0.5 cm FWHM [31, 33], but no depth resolution was given for 
the 2 cm cubic crystal. In addition, the geometry of crystals 
will determine their spatial resolution, the specific parameters 
of which can be seen in Sec. IIB of the paper. The organic 
plastic scintillator ACS used an energy resolution of 10 keV 
(lo Gaussian), a trigger threshold of 100 keV, and a detection 
efficiency of 99.9% for charged particles. 


Since the MEGAIib package provides the function of user- 
defined particle sources, we can flexibly set the particle 
sources for the simulation according to our needs. For our 
purposes, we have used monochromatic, negative power law, 
and the file format (generated by the LEOBackground soft- 
ware package) for the energy spectrum, as well as beam 
parameters for the far-field point source (i.e., homogeneous 
beam) and the far-field area source. When setting up parti- 
cle sources, we also took into account the fact that far-field 
sources are so far away and that they arrive at the detector in 
the form of plane waves. In our simulation experiments, a 
total of three types of sources were applied to the three sim- 
ulations. First, a monochromatic homogeneous beam with 
a flux of 1.0 ph/cm?/s was simulated and used to irradiate 
the mass model of the polarimeter. The polarization response 
of the polarimeter was studied by varying the energy, an- 
gle, and polarization direction of incident photons. Next, 
the homogeneous beam was repeated with a power law en- 
ergy spectrum to simulate a discrete celestial source of lin- 
early polarized gamma rays (i.e., a Crab-like source). For the 
Crab-like source, events were generated by linearly polarized 
and unpolarized photon beams with an energy spectrum of 
4x1073 E72 ph/cm?/s/MeV between 0.1 and 10 MeV [45]. 
A simulation experiment was performed using this type of 
source to investigate the polarization response of the po- 
larimeter to realistically polarized photons and its polariza- 
tion performance. Third, an isotropic beam source (i.e., a 
far-field area source) with a spectrum in file form was con- 
structed. In fact, this type of source was used to simulate the 
space orbit environment. Since the polarimeter is expected 
to complete a satellite mission, it is quite necessary to con- 
sider the impact of the complex space orbit environment on 
the detector. The background environment depends to a large 
extent on the orbit in which the satellite is operating. In this 
case, a typical LEO with an altitude of 550 km and a 0° incli- 
nation was selected. The detection of the corresponding or- 
bital background environments by the Beppo-SAX [46] and 
AGILE [47] missions has made orbital background environ- 
ments like this one well-known. Figure 6 shows the back- 
ground environment used in the polarimeter simulation ex- 
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periments, calculated and given by LEOBackground. As can 
be seen in Fig. 6, the background components consist of cos- 
mic and albedo photons, hadrons (e.g., neutrons, protons, and 
alphas), and leptons (e.g., electrons and positrons). It is worth 
mentioning that the energy spectra in file form that we used 
for the simulation experiments are exactly these background 
energy spectra. 
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Fig. 6. (Color online) The background environment of the polarime- 
ter on an orbit with an altitude of 550 km and a 0° inclination. 


For data analysis, first, a coincident event filter was adopted 
to select events that satisfy at least one hit in the converter 
as well as at least one hit in the absorber. Next, the Comp- 
ton sequence reconstruction (CSR) algorithm with trajectory 
tracking [37] was applied to these selected events for Comp- 
ton event reconstruction, which would require that the first hit 
in the event needs to occur in the converter. Then the infor- 
mation related to the incident photon could be obtained, such 
as the initial energy, the incident direction, and the scatter- 
ing angle of the incident photon. All of the above was done 
by Revan contained in MEGAIib, where the event selection 
method ensures the integrity of Compton events and plays a 
vital role in reducing background events. Finally, the recon- 
structed events would be further processed and analyzed by 
the Mimrec tool in MEGAIib, including event cuts, plotting of 
the reconstructed energy spectrum, the evaluation of the an- 
gular resolution measurement (ARM) specified as the angular 
distance between a known source position and the closest re- 
constructed position on the Compton cone, calculation and 
correction of modulation curves, image reconstruction, etc. 
When we analyzed the reconstructed events, we adopted the 
same and fixed event cuts for events generated by monoen- 
ergetic photons, including a +30 photopeak energy window 
and a +30 ARM cut. In addition, there was a variable event 
cut, namely the scattering angle cut, which would be adjusted 
according to the different photon energies, and the scatter- 
ing angle window is roughly chosen between 40° and 110°. 
For events generated by Crab-like sources, both a +30 ARM 
window and a [60°, 110°] Compton scattering angle window 
were applied as event cuts, and no energy cut was used. Sim- 
ilarly, for events generated by the background sources in the 
polarimeter, the same event analysis method was used. After 
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the above data processing series, the distribution of azimuthal 
scattering angles was obtained. Then, the polarization infor- 
mation of the incident photons would be gained from the pa- 
rameters fitted to the azimuthal distribution. 


Ill. RESULTS AND DISCUSSION 


In simulated experiments, various radiation physics pro- 
cesses occur when soft gamma rays are incident in the po- 
larimeter. These physical processes are mainly photoelec- 
tric effect, Compton scattering (see yı and y2 in Fig. 7), 
and electron-positron pair creation, with Compton scattering 
dominating. Of course, gamma rays will also not interact with 
matter at all, such as 73, ya, and 75 in Fig. 7. As the polarime- 
ter designed in this paper is based on the Compton scattering 
principle for the detection of linearly polarized gamma rays, 
the Compton scattering process is of interest to us. We re- 
quire that the incident gamma rays can satisfy the condition 
that they first undergo Compton scattering in the converter 
and are finally stopped by the absorber. The typical effec- 
tive Compton scattering process and the detection process of 
the polarimeter are described as follows. Generally, when a 
linearly polarized photon is incident in the polarimeter, the 
photon first undergoes Compton scattering in the converter. 
At this time, the information on the position of the scattering 
point and the energy of the recoil electron can be measured 
by the converter, while the recoil electron will be tracked if 
it creates a trail in the converter. Then the absorber absorbs 
the scattered photon and obtains the energy of the scattered 
photon and the position information of the absorption point. 
Of course, in addition to the above, there are also cases where 
recoil electrons pass through the converter and are absorbed 
by the absorber (see y2 in Fig. 7), as well as cases where scat- 
tered photons undergo multiple Compton scattering before 
being completely stopped by the absorber (see ys in Fig. 7). It 
is clear from Fig. 7 that yı and 2 are valid events that satisfy 
the condition. Finally, the simulation results can be obtained 
by analyzing these valid events. 


A. Polarization response to monoenergetic photons 


The fully linearly polarized (100%) and non-polarized 
(0%) homogeneous on-axis photon beams with an energy of 
200 keV were used in the simulation experiments, where the 
other configuration parameters of both beams were identical 
(after here, photon beams are considered to be on-axis inci- 
dence if not otherwise specified). For a 100% linearly polar- 
ized photon beam, the polarization vector was (1, 0, 0). The 
raw azimuthal scatter angle distributions were obtained after 
the analysis of simulated data, as shown in Fig. 8. As can 
be seen in Fig. 8, the polarimeter has a significant response 
to polarized photons (see the blue curve) relative to unpo- 
larized photons (see the black curve). It should be noted, 
however, that although the blue curve approximately obeys 
the cosine distribution, it is affected and distorted by the sys- 
tematic modulation due to the non-uniformity of the detector, 
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Fig. 7. (Color online) Interaction processes and traces left by gamma 
rays with an energy of 1000 keV incident on-axis into the polarime- 
ter. Note: To make the trajectory of the particles in the polarimeter 
clearer, we have given a front view of the detector. 


and it might be difficult to be fitted by Eq. (3). In addition, 
for the azimuthal angle distribution generated by unpolarized 
photons in the detector (see the black curve in Fig. 8), an ap- 
parent spurious modulation (the black curve is not perfectly 
uniformly distributed) is shown owing to the non-uniformity 
of the detector response along the azimuthal angle (i.e., the 
effect of the non-uniform response). Therefore, the raw az- 
imuthal angle distribution needs to be corrected in order to 
eliminate the impact of this systematic modulation on the 
modulation curve. To obtain the corrected azimuthal angle 
distribution, we can follow the formula [2, 37] 


f(y) 


fron (p) 


to correct the raw distribution. Here, f(y) is the raw az- 
imuthal angle distribution generated by polarized photons (as 
the blue curve in Fig. 8), fron (p) is the azimuthal angle dis- 
tribution produced by unpolarized photons with the same en- 
ergy and incident direction as the polarized photon (as the 
black curve in Fig. 8), and f.o,(v) is the corrected azimuthal 
angle distribution. Figure 9 illustrates the corrected modula- 
tion curve (after here, all modulation curves shown and an- 
alyzed are corrected). Compared with the raw modulation 
curve (see the blue curve in Fig. 8), the corrected modula- 
tion curve obeys the cosine distribution perfectly and can be 
well-fitted by Eq. (3) (see Fig. 9). The modulation factor 
H100 can be calculated using the fitting parameters, which is 
0.80+0.01. Also, 0.2°+0.2° is obtained to serve as the po- 
larization direction of the incident photons, which is in good 
agreement with the actual direction (0°). From the above sim- 
ulation results, it can be evident that the polarization response 
and performance of the polarimeter are excellent. 

Figure 10 presents the corrected modulation curves fitted 
by Eq. (3) for four different polarization directions of the pho- 
ton beam with an energy of 200 keV. The polarization pa- 
rameters, such as polarization angle and modulation factor, 
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Fig. 8. (Color online) The raw azimuthal scatter angle distributions 
obtained from simulations that 100% (blue) and 0% (black) linearly 
polarized photon beams with an energy of 200 keV were on-axis 
incidence. 


6000 


5000 


4000 


3000 


Corrected counts 


2000 


1000 


Pireieletaiterrr tives lair derevbiviales at 
430 -135 -90 —-45 0 45 90 135 180 
Azimuthal scatter angle (degree) 


Fig. 9. (Color online) The azimuthal scatter angle distribution was 
obtained after the correction of the raw distribution (blue) in Fig. 8. 
The red curve was generated by fitting the corrected distribution by 


Eq. (3). 


of the incident photons with four different polarization direc- 
tions were obtained and listed in Table 2. From the simulation 
results in Table 2, it can be concluded that the polarimeter 
can precisely determine the polarization angle with an error 
of 0.2° no matter how the polarization direction of the inci- 
dent photons is changed, and the obtained modulation factors 
are quite large and remain essentially constant. 

During the data analysis, the Compton scattering events 
produced by polarized photons with an energy of 200 keV 
were selected with the Compton scattering angle as the filter- 
ing condition, and the dependence of the modulation factor 
on the Compton scattering angle can be obtained as shown in 
Fig. 11(a). As can be seen in Fig. | 1(a), the modulation factor 
is smaller in the region of large and small scattering angles, 
yet it increases significantly in the region of medium scatter- 
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Fig. 10. (Color online) Corrected modulation curves and their best- 
fit curves (red). The polarization directions of the photon beams are 
0° (black), -60° (green), -90° (blue), and -180° (magenta), respec- 
tively. 


TABLE 2. Polarization parameters for polarized photons with polar- 
ization directions 0°, -60°, -90°, and -180°, respectively. 


Polarization direction Polarization angle Modulation factor 


setup simulated H100 
(degree) (degree) simulated 
0 0.20.2 0.80+0.01 
-60 -60.2+0.2 0.80+0.01 
-90 -90.2+0.2 0.80+0.01 
-180 -179.8+0.2 0.80+0.01 


ing angles. Figure | 1(b) shows the relationship between the 
energy of the incident polarized photons (100 keV, 200 keV, 
300 keV, 500 keV, 800 keV, 1000 keV, 2000 keV, 3000 keV, 
5000 keV, and 10000 keV) and the modulation factor. It is ap- 
parent that the modulation factor decreases with the increase 
of the incident photon energy. From the simulation results in 
Fig. 11, we can conclude that the modulation of the azimuthal 
angle distribution is most significant at lower energies and 
medium Compton scattering angles, which is consistent with 
the theoretical results in Fig. 3. 


Polarized photons with energies of 200 keV, 500 keV, and 
1000 keV were simulated to be incident into the polarimeter 
at five different polar angles (0°, 20°, 40°, 60°, and 80°). Then 
three curves of the modulation factor versus incident polar an- 
gle were obtained, as shown in Fig. 12. The three relationship 
curves show that the modulation factor is almost unaffected 
by the off-axis incidence (not dependent on the incident po- 
lar angle of the photon), indicating that the polarimeter can 
still maintain a stable and remarkable polarization response 
to photons with different incidence angles. In addition, the 
results verify that such a polarimeter design has a large FoV. 
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Fig. 11. (Color online) Dependence of the modulation factor on 
the Compton scattering angle (left panel, (a)) and the photon en- 
ergy (right panel, (b)) for the monoenergetic photons. 


1= 
0.8— e—a 
a F 
= pe — + — l 
g o6 
8 L 
sab eee ee 
= 04+ 
3 L 
© 
£ B 
0 2L —e— Energy: 200keV 
al —a Energy: 500keV 
L —— Energy: 1000keV 


tirehrerrrlrrrihae rritin ds 
$ 0 10 20 30 40 50 60 70 80 90 


y-Ray incident Polar angle (degree) 


Fig. 12. (Color online) Dependence of the modulation factor on the 
incident polar angle for photons with energies of 200 keV, 500 keV, 
and 1000 keV. 


B. Polarization response to a Crab-like source 


Two Crab-like sources with a fully linearly polarized 
(100%) gamma-ray beam and an unpolarized (0%) gamma- 
ray beam (see Sec. IID for detailed parameters) were used 
in the simulation experiments of on-axis incidence on the 
polarimeter, provided that all other configurations were the 
same. The modulation factor was obtained by fitting the cor- 
rected modulation curve generated by 100% linearly polar- 
ized photons with Eq. (3), which had a value of 0.76+0.01. 
It can be observed that the polarimeter has a good response 
to continuous polarization photon spectra in addition to being 
sensitive to monoenergetic polarization photon spectra. 

Since the polarization degree of the polarized photons 
emitted by the cosmic gamma-ray sources may not always be 
100% but may be partially polarized, it is necessary to study 
the polarization response of the polarimeter to photons of dif- 
ferent polarization degrees. The simulated modulation factor 
Lt (see Eq. (4)) and the simulated polarization degree P (see 
Eq. (5)) were obtained by illuminating the polarimeter with 
photon beams emitted by the Crab-like sources with differ- 
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ent polarization degrees (1.0, 0.8, 0.6, 0.4, 0.2, and 0.1). The 
simulation results for the two polarization parameters men- 
tioned above are presented in Table 3. From the data in the 
first and third columns in Table 3, it can be found that the 
modulation factor decreases as the polarization degree of the 
photon falls, which indicates that the sensitivity of the de- 
tector to photons with small polarization degree decreases, 
and this phenomenon is consistent with the theory. In addi- 
tion, the data in the first and second columns in Table 3 show 
that the simulated polarization degrees are in good agreement 
with the actual set polarization degrees within the error tol- 
erance. Therefore, we can assume that the polarimeter has a 
good ability to preserve the polarization degree of the inci- 
dent photons or that the polarimeter can recover the incident 
photon polarization degree well. 


TABLE 3. The simulated modulation factor u (see Eq. (4)) and the 
simulated polarization degree P (see Eq. (5)) for the different polar- 
ization degrees of the incident photons from the Crab-like sources. 


Crab-like source Polarization degree Modulation factor 


polarization degree P H 

setup simulated simulated 

1 1.00+0.02 0.76+0.01 

0.8 0.80+0.02 0.61+0.01 

0.6 0.59+0.02 0.45+0.01 

0.4 0.39+0.01 0.30+0.01 

0.2 0.20+0.01 0.15+0.01 

0.1 0.09+0.01 0.07+0.01 


The polarization sensitivity (MDP) can be defined as 
Eq. (6), a crucial performance parameter for a polarimeter, 
and characterizes the minimum polarization that a polarime- 
ter can detect. The MDP of the polarimeter can be estimated 
by simulation. When the Crab-like source and background 
components were used for the simulation experiments of the 
polarimeter, the relationship curve between the MDP@30 
and the observation time was obtained, as shown in Fig. 13. 
It can be seen that the polarization detection capability of 
the polarimeter increases with observation time. During the 
twenty-four hours observation period, the MDP@3c of the 
polarimeter is 8.2%. When the observation time reaches 10° 
seconds (~277.8 hours), the polarimeter is capable of detect- 
ing a 2.4% polarized Crab-like source. In conclusion, the po- 
larimeter with excellent polarization detection capability can 
be well used to detect the polarization information of soft 
gamma rays in space to study various astronomical sources. 


IV. SUMMARY 


As a new astronomical tool, cosmic gamma-ray polar- 
ization can help explain various astrophysical phenomena 
and mechanisms from a unique dimension, yet it isn’t easy 
to measure. Currently, there is no effective dedicated soft 
gamma-ray polarimeter operating in orbit. Small satellite 
technology is developing rapidly and has unique advantages 
compared to traditional large satellites. Given those men- 
tioned above, developing a high-performance dedicated space 
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Fig. 13. (Color online) The polarization sensitivity (MDP) of the 
polarimeter for a Crab-like source as a function of the observation 
time. 


soft gamma-ray polarimeter is necessary. The idea of using 
a constellation of multiple microsatellites to detect the polar- 
ization of gamma rays in space is promising. In this paper, 
we have studied one of the microsatellites in the constella- 
tion in detail. First, a novel structure was designed, which in- 
cluded a multilayer double-sided silicon strip converter in the 
upper part of the detector, a pixel-type CsI array absorber on 
five sides except the upper side of the converter, and a plastic 
scintillator ACS at the outermost side. Then, the mass model 
of the polarimeter was constructed and simulated in detail 
by setting the performance parameters of each sub-detector, 
the different gamma-ray polarization sources, and the back- 
ground components using the dedicated simulation software 
based on the Monte Carlo method. For 100% polarized pho- 
tons with an energy of 200 keV that were on-axis incidence, 
the polarimeter achieved a modulation factor of 0.80+0.01. 
The modulation factor of the polarimeter showed little de- 
pendence on the direction of incidence of the photons. In 
addition, the polarimeter could precisely measure the polar- 
ization angle of the incident photon at an accuracy of 0.2 de- 
grees. For on-axis incident photons in the energy range of 
0.1-10 MeV from a Crab-like polarization source, the polar- 
ization response of the polarimeter remained remarkable with 
a modulation factor of 0.760.01 (100% polarization), and 
the polarization degree of the incident photons could be accu- 
rately measured with a maximum error of 0.02. At the same 
time, the polarimeter exhibited excellent polarization sensi- 
tivity, and the MDP @3c of the polarimeter could reach 2.4% 
at an exposure of 10° seconds. Undeniably, the outstanding 
performance demonstrated by the polarimeter also depends to 
some extent on its wide FoV (~27 sr) and potent background 
suppression. Additionally, the simulation results of the po- 
larimeter that the polarization response of the polarimeter is 
most pronounced at lower energy and medium scattering an- 
gle were in good agreement with the theoretical law. 


In summary, the idea of using microsatellites to detect soft 


717 gamma-ray polarization proposed in this paper is promising 
zs and can bring new thoughts for future space exploration. Ad- 
719 ditionally, the polarimeter of the novel configuration designed 
72 in this paper has shown excellent polarization detection capa- 
721 bility after being verified by simulation experiments, which 
72 has great potential to be competent for the future task of soft 
73 gamma-ray polarization detection in space. In this paper, the 


[1] R. Mignani, A. Shearer, A. Slowikowska et al., Astronomi- 
cal Polarisation from the Infrared to Gamma Rays, (Springer, 
Switzerland Cham, 2019), pp. 109—111. 

[2] S. Komura, A. Takada, Y. Mizumura et al., Imaging Po- 
larimeter for a Sub-MeV Gamma-Ray All-sky Survey Using 
an Electron-tracking Compton Camera. Astrophys. J. 839, 41 
(2017). doi: 10.3847/1538-4357/aa68dc 

[3] K. Toma, T. Sakamoto, B. Zhang et al., Statistical Properties 

of Gamma-Ray Burst Polarization. Astrophys. J. 698, 1042— 

1053 (2009). doi: 10.1088/0004-637X/698/2/1042 

V. Tatischeff, A.D. Angelis, C. Gouiffés et al, e- 

ASTROGAM mission: a major step forward for gamma-ray 

polarimetry. J. Astron. Telesc. Inst. 4, 011003 (2017). doi: 

10.1117/1.JATIS.4.1.011003 

N. Produit, F. Barao, S. Deluit et al., POLAR, a compact detec- 

tor for Gamma Ray Bursts photon polarization measurements. 

Nucl. Instrum. Methods Phys. Res. Sect. A 550, 616—625 

(2005). doi: 10.1016/j.nima.2005.05.066 

Y.Z. Fan, Interpretation and implications of the non-detection 

of GeV spectrum excess by the Fermi Gamma-ray Space 

Telescope in most gamma-ray bursts. Mon. Not. R. As- 

tron. Soc. 397, 1539—1548 (2009). doi: 10.1111/j.1365- 

2966.2009.15018.x 

A.D. Angelis, V. Tatischeff, I.A. Grenier et al., Science with 

e-ASTROGAM: A space mission for MeV—GeV gamma-ray 

astrophysics. J. High Energy Astrop. 19, 1—106 (2018). doi: 
10.1016/j.jheap.2018.07.001 

J.D. Schnittman, J-H. Krolik, X-ray Polarization from Accret- 

ing Black Holes: Coronal Emission. Astrophys. J. 712, 908— 

924 (2010). doi: 10.1088/0004-637X/7 12/2/908 

[9] F. Lei, A.J. Dean, G.L. Hills, Compton Polarimetry in Gamma- 
Ray Astronomy. Space Sci. Rev. 82, 309—388 (1997). doi: 
10.1023/A:1005027107614 

[10] H.C. Zhang, C. Diltz, M. Böttcher, Radiation and Polarization 

Signatures of 3D Multi-zone Time-dependent Hadronic Blazar 
Model. Astrophys. J. 829, 69 (2016). doi: 10.3847/0004- 
637X/829/2/69 

[11] M. Böttcher, Progress in Multi-Wavelength and Multi- 

769 Messenger Observations of Blazars and Theoretical Chal- 

770 lenges. Galaxies 7, 20 (2019). doi: 10.3390/galaxies7010020 

77 [12] P. Laurent, J. Rodriguez, J. Wilms et al., Polarized Gamma-Ray 

Emission from the Galactic Black Hole Cygnus X-1. Science 

773 332, 438—439 (2011). doi: 10.1126/science. 1200848 

774 [13] H. Krawczynski, A. Garson, Q. Guo et al., Scientific prospects 

for hard X-ray polarimetry. Astropart. Phys. 34, 550—567 

776 (2011). doi: 10.1016/j.astropartphys.2010.12.001 

777 [14] M.C. Weisskopf, E.H. Silver, H.L. Kestenbaum et al., A preci- 

778 sion measurement of the X-ray polarization of the Crab Nebula 

779 without pulsar contamination. Astrophys. J. 220, L117—L121 

780 (1978). doi: 10.1086/182648 

731 [15] M.L. McConnell, J.M. Ryan, D.M. Smith et al., RHESSI as 

782 a Hard X-Ray Polarimeter. Sol. Phys. 210, 125—142 (2002). 


[4] 


[5 


=a 


[6 


= 


[7] 


[8] 


760 
761 
762 
763 
764 
765 
766 
767 
768 


772 


775 


11 


724 configuration design, simulation methods, and simulation re- 
s sults of the polarimeter will lay a solid foundation for the de- 
velopment of the future polarization detection satellite proto- 
7 type, will strongly promote the smooth development of possi- 
zs ble future constellation programs, and will provide a valuable 
729 reference for the design of other polarization detectors in the 


730 future. 


7: 


N 


726 


7 


N 


doi: 10.1023/A:1022413708738 

M. Forot, P. Laurent, I. Grenier et al., Polarization of the 
Crab Pulsar and Nebula as Observed by the INTEGRAL/IBIS 
Telescope. Astrophys. J. Lett. 688, L29—L32 (2008). doi: 
10.1086/593974 

D. Yonetoku, T. Murakami, S. Gunji, et al., DETECTION OF 
GAMMA-RAY POLARIZATION IN PROMPT EMISSION 
OF GRB 100826A. Astrophys. J. Lett. 743, L30 (2011). doi: 
10.1088/2041-8205/743/2/L30 

H.L. Xiao, Y.W. Dong, B.B. Wu et al., POLAR Gamma Ray 
Burst Polarimeter onboard TG-2 Spacelab. Manned Spaceflight 
21, 32-36+43 (2015). (in Chinese) doi: 10.3969/j.issn.1674- 
5825.2015.01.006 

V. Bhalerao, D. Bhattacharya, A. Vibhute et al., The Cadmium 
Zinc Telluride Imager on AstroSat. J. Astrophys. Astron. 38, 
31 (2017). doi: 10.1007/s12036-017-9447-8 

Y. Kishimoto, S. Gunji, Y. Ishigaki et al., Basic Per- 
formance of PHENEX: A Polarimeter for High ENErgy 
X rays. IEEE T. NUCL. SCI. 54, 561-566 (2007). doi: 
10.1109/TNS.2007.897827 

S. Larsson, M. Pearce, PoGO : The polarised gamma-ray ob- 
server. Nucl. Instrum. Methods Phys. Res. Sect. A 525, 148- 
152 (2004). doi: 10.1016/j.nima.2004.03.036 

M. Friis, M. Kiss, V. Mikhalev et al., The PoGO+ Balloon- 
Borne Hard X-ray Polarimetry Mission. Galaxies 6, 30 (2018). 
doi: 10.3390/galaxies6010030 

P.F. Bloser, J.S. Legere, J.R. Macri et al., GRAPE - A Balloon- 
Borne Gamma-Ray Polarimeter Experiment. Chin. J. Astron. 
Astrophys. 6, 393—397 (2006). doi: 10.48550/arXiv.astro- 
ph/05083 14 

T. Satoru, A. Shigeki, K. Keiki et al., GRAINE project: 
The first balloon-borne, emulsion gamma-ray telescope ex- 
periment. Prog. Theor. Exp. Phys. 2015, 43H01 (2015). doi: 
10.1093/ptep/ptv046 

P.F. Bloser, J.M. Ryan, M.L. McConnell et al., The 
MEGA project: Science goals and hardware devel- 
opment. New Astron. Rev. 50, 619-623 (2006). doi: 
10.1016/j.newar.2006.06.001 

J.-L. Chiu, S.E. Boggs, H.-K. Chang et al., The upcoming 
balloon campaign of the Compton Spectrometer and Imager 
(COSI). Nucl. Instrum. Methods Phys. Res. Sect. A 784, 359- 
363 (2015). doi: 10.1016/j.nima.2014.11.099 

X. Wen, X.W. Zhang, Analysis of current developments in 
the field of modern small satellites. Journal of Shenyang 
Aerospace University 30, 33-36+46 (2013). (in Chinese) doi: 
10.3969/j.issn.2095-1248.2013.05.007 

Z.G. Bai, Development Achievements and Prospects of 
China Modern Small Satellite. Spacecraft Engineering 
28, 1-8 (2019). (in Chinese) doi:  10.3969/j.issn.1673- 
8748.2019.02.001 

J.X. Wen, X.T. Zheng, J.D. Yu et al., Compact CubeSat 
Gamma-ray detector for GRID mission. Nucl. Sci. Tech. 32, 


783 
784 [16] 
785 
786 
787 
788 
789 
790 


[17] 


791 
792 [18] 
793 
794 
795 
796 
797 
798 
799 
800 


[19] 


[20] 


801 
802 
803 
804 


[21] 


805 
soe [22] 
807 
808 
809 
810 


[23] 


811 
812 
813 
814 


[24] 


815 
816 
817 
818 
819 
820 
821 
822 
823 
824 


[25] 


[26] 


825 [27] 
826 
827 
828 
829 
830 
831 
832 
833 
834 


[28] 


[29] 


m 857 


835 107-117 (2021). doi: 10.1007/S41365-021-00937-4 
s36 [30] O. Engvold, J.-C. Vial, A. Skumanich, The Sun as a Guide to 
837 Stellar Physics , (Elsevier, Amsterdam, 2019), pp. 189. 
838 [31] W.F. Liu, X.M. Liu, S.W. Tang et al., Study of the detection 
839 unit of the calorimeter for the next generation Compton tele- 
840 scope. Nuclear Techniques 43, 20-28 (2020). (in Chinese) doi: 
10.11889/j.0253-3219.2020.hjs.43.010203 
V. Bindi, A.D. Guerra, G. Levi et al., Preliminary study of sil- 
icon photomultipliers for space missions. Methods Phys. Res. 
Sect. A 572, 662—667 (2007). doi: 10.1016/j.nima.2006.12.011 
X.M. Liu, S.W. Tang, W.F. Liu et al., Study of the y -Ray De- 
tector of Dual-ended Readout Based on SiPM. Nuclear Physics 
Review 37, 757—764 (2020). (in Chinese) doi: 10.11804/Nu- 
clPhysRev.37.2019CNPC03 
A. Zoglauer, R. Andritschke, F. Schopper, MEGAILib — The 
Medium Energy Gamma-ray Astronomy Library. New Astron. 
Rev. 50, 629—632 (2006).doi: 10.1016/j.newar.2006.06.049 
[35] R. Brun, F. Rademakers, ROOT — An object oriented data 
853 analysis framework. Nucl. Instrum. Methods Phys. Res. Sect. 
854 A 389, 81—86 (1997). doi: 10.1016/S0168-9002(97)00048-X 
855 [36] S. Agostinelli, J. Allison, K. Amako et al., Geant4—a simu- 
lation toolkit. Nucl. Instrum. Methods Phys. Res. Sect. A 506, 
250-303 (2003). doi: 10.1016/S0168-9002(03)01368-8 
68 [37] A.C. Zoglauer, First Light for the Next Generation of Comp- 
ton and Pair Telescopes. Dissertation for the Doctoral Degree 
(Technical University of Munich, Munich, 2005), pp. 41—125. 
861 [38] R. Rando, The all-sky medium energy gamma-ray obser- 
vatory. J. Instrum. 12, C11024 (2017). doi: 10.1088/1748- 
» 863 0221/12/11/C11024 
= s64 [39] V. Schénfelder, The imaging gamma-ray telescope COMP- 


841 
842 [32] 
843 
844 
s45 [33] 
846 
847 
848 
s49 [34] 
850 

851 


852 


856 


859 
860 


862 


865 
866 
867 
868 
869 
870 


[40] 


[41] 
871 
872 
873 
874 
875 


[42] 


876 
877 
878 
879 
880 


[43] 


881 
882 
883 
884 
885 [45] 
886 
887 
sess [46] 
889 
890 
s91 [47] 
892 

893 


[44] 


12 


TEL aboard GRO. Adv. Space Res. 11, 313-322 (1991). doi: 
10.1016/0273-1177(91)90183-K 

S.E. Boggs, The Advanced Compton Telescope mis- 
sion. New Astron. Rev. 50, 604-607 (2006) doi: 
10.1016/j.newar.2006.06.076 

K. Kamiya, Compton Recoil Electron Tracking With the TIGRE 
Gamma-Ray Balloon Experiment. Dissertation for the Doctoral 
Degree (University of California, Riverside, Riverside, 2011), 
pp. 53—61. 

M. Galloway, A. Zoglauer, S.E. Boggs et al, A com- 
bined Compton and coded-aperture telescope for medium- 
energy gamma-ray astrophysics. A&A 614, A93 (2018) doi: 
10.1051/0004-6361/201731122 

A. Zoglauer, S.E. Boggs, M. Galloway et al., Design, im- 
plementation, and optimization of MEGAIib’s image re- 
construction tool Mimrec. A&A 652, 568-571 (2011) doi: 
10.1016/j.nima.2010.08.043 

P. Cumani, M. Hernanz, J. Kiener et al., Background for a 
gamma-ray satellite on a low-Earth orbit. Exp. Astron. 47, 273- 
—302 (2019). doi: 10.1007/s10686-019-09624-0 

A. Akyüz, D. Bhattacharya, D.D. Dixon et al., Polarization 
measurements with an imaging gamma-ray telescope. Exp. As- 
tron. 6, 275—284 (1995). doi: 10.1007/BF00418325 

R. Campana, M. Orlandini, E.D. Monte et al., The radiation 
environment in a low earth orbit:the case of BeppoSAX. Exp. 
Astron. 37, 599-613 (2014). doi: 10.1007/s10686-014-9394-1 
M. Tavani, G. Barbiellini, A. Argan et al., The AGILE 
mission. Astron. Astrophys. 502, 995—1013 (2009). doi: 
10.105 1/0004-636 1/2008 10527 


